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We present a simple quantum network, based on the controUed-SWAP gate, that can extract 
certain properties of quantum states without recourse to quantum tomography. It can be used used 
as a basic building block for direct quantum estimations of both linear and non-linear functionals 
of any density operator. The network has many potential applications ranging from purity tests 
and eigenvalue estimations to direct characterization of some properties of quantum channels. Ex- 
perimental realizations of the proposed network are within the reach of quantum technology that is 
currently being developed. 
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Certain properties of a quantum state g, such as its 
purity, degree of entanglement, or its spectrum, are of 
significant importance in quantum information science. 
They can be quantified in terms of linear or non-Unear 
functionals of g. Linear functionals, such as average val- 
ues of observables {A}, given by TrAgi, are quite com- 
mon as they correspond to directly measurable quanti- 
ties. Non-linear functionals of state, such as the von 
Neumann entropy —"Tiglng, eigenvalues, or a measure 
of purity Trp^, are usually extracted from g by classical 
means i.e. g is first estimated and once a sufficiently pre- 
cise classical description of g is available, classical evalu- 
ations of the required functionals can be made. However, 
if only a limited supply of physical objects in state g is 
available, then a direct estimation of a specific quantity 
may be both more efficient and more desirable [Q. For 
example, the estimation of purity of g does not require 
knowledge of all matrix elements of g, thus any prior 
state estimation procedure followed by classical calcula- 
tions is, in this case, inefficient. However, in order to 
bypass tomography and to estimate non-linear function- 
als of g more directly, we need quantum networks |^ 
performing quantum computations that supersede clas- 
sical evaluations. 

In this paper, we present a simple quantum network 
which can be used as a basic building block for direct 
quantum estimations of both linear and non-linear func- 
tionals of any g. The network can be realized as mul- 
tiparticle interferometry. While conventional quantum 
measurements, represented as quantum networks or oth- 
erwise, allow the estimation of Ag for some Hermitian 
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FIG. 1: Both the visibility and the shift of the interfer- 
ence patterns of a single qubit (top line) are affected by the 
controlled-(7 operation on a general state, p. 



A, our network can also provide a direct estimation of 
the overlap of any two unknown quantum states ga and 
gb, i.e. Tt: gagb- Here, and in the following, we use ter- 
minology developed in quantum information science. For 
a comprehensive overview of this terminology, including 
quantum logic gates and quantum networks see, for ex- 
ample, [§. 

In order to explain how the network works, let us start 
with a general observation related to modifications of 
visibility in interferometry. Consider a typical interfero- 
metric set-up for a single qubit: Hadamard gate, phase 
shift (p, Hadamard gate, followed by a measurement in 
the computational basis. We modify the interferome- 
ter by inserting a controUed-C/ operation between the 
Hadamard gates, with its control on the qubit and with U 
acting on a quantum system described by some unknown 
density operator p, as shown in Fig. ^ The action of 
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FIG. 2: A quantum network for direct estimations of both 
linear and non-linear functions of state. The probability of 
finding the control (top line) qubit in state | 0) at the output 
depends on the overlap of the two target states (two bottom 
lines). Thus estimation of this probability leads directly to 
an estimation oiTx QaQb = v ~ 2Po — 1. 



In general, in order to evaluate the spectrum of any 
d X d density matrix g we need to estimate d — 1 pa- 
rameters Tig^, Tig^,... Trg^. For this we need the 
controUcd-shift operation, which is a generalization of the 
controUed-swap gate. Given k systems of dimension d we 
define the shift as 
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(3) 



for any pure states | 0) . Such an operation can be easily 
constructed by cascading fc — 1 swaps V. If we extend 
the network and prepare p — g'^'' at the input then the 
interference will be modified by the visibility factor. 
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the controUed-C/ on p modifies the interference pattern 
by the factor ||] , 
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where v is the new visibility and a is the shift of the 
interference fringes, also known as the Pancharatnam 
phase Q. The observed modification of the visibility 
gives an estimate of Tr Up, i.e. the average value of the 
unitary operator U in state p. Let us mention in passing 
that this property, among other things, allows to esti- 
mate an unknown quantum state p as long as we can 
estimate TvUkP for a set of unitary operators Uk which 
form a basis in the vector space of density operators. 

Let us now consider a quantum state p of two separa- 
ble subsystems, such that p = ga® gb- We choose our 
controUed-C/ to be the controlled-^, where V is the swap 
operator, defined as, V \ a)^ | (3) g = | | a)^, for any 
pure states |q;)^ and 1/3)^. In this case, the modifica- 
tion of the interference pattern given by Eq. (|]) can be 
written as, 



V = Tr {ga ® Qb) = Tr gagb- 



(2) 



which is easily proved using the spectral decomposition 
of ga and gi,. Since Tr gagb is real, we can fix (ys = and 
the probability of finding the qubit in state | 0) at the 
output, Pq, is related to the visibility by v = 2Po — 1- 
This construction, shown in Fig. (^), provides a direct 
way to measure Tr gagb (c.f. for a related idea). 

There are many possible ways of utilizing this result. 
For pure states ga ~ \a) {a\ and gb — \ f3) {f3\ the formula 
above gives Tr gagb — | ( a | /3 ) ^ i.e. a direct measure of 
orthogonality of | a) and \ f3). If we put ga = Qb = Q then 
we obtain an estimation of Tr g^ . In the single qubit case, 
this measurement allows us to estimate the length of the 
Bloch vector, leaving its direction completely undeter- 
mined. For qubits Tr g^ gives the sum of squares of the 
two eigenvalues which allows to estimate the spectrum of 
g. 



Thus measuring the average values of V^^'' for k = 2, 3...d 
allows us to evaluate the spectrum of ^) Although 
we have not eliminated classical evaluations, we have re- 
duced them by a significant amount. The average values 
of for k = 2,3...d provide just enough information 
to evaluate the spectrum of g but certainly not enough 
to estimate the whole density matrix. 

So far we have treated the two inputs, ga and gb in 
a symmetric way. However, there are several interest- 
ing applications in which one of the inputs, say go, is 
predetermined and the other is unknown. For example, 
projections on a prescribed vector \ '4'), or on the sub- 
space perpendicular to it, can be implemented by choos- 
ing ga — \ {'^\- By changing the input state | ijj) we 
effectively "reprogram" the action of the network which 
then performs different projections. This property can 
be used in quantum communication, in a scenario where 
one carrier of information, in state | ip) , determines the 
type of detection measurement performed on the second 
carrier. N.B. as the state \tp) oia single carrier cannot be 
determined, the information about the type of the mea- 
surement to be performed by the detector remains secret 
until the moment of detection. 

Another interesting application is the estimation of the 
extremal eigenvalues and eigenvectors of gb without re- 
constructing the entire spectrum. In this case, the input 
states are of the form | V') (V' l®P6 f^nd we vary | tjj) search- 
ing for the minimum and the maximum of v = ('0 I ^?b I "0) ■ 
This, at first sight, seems to be a complicated task as it 
involves scanning 2(c? — 1) parameters of ip. The visibility 
is related to the overlap of the reference state | 0) and gb 
by, 

vv- = Tr (^|V^)(V'|^A,|77z)(r/.|^ 

i i 

where J^iPi — 1- This is a convex sum of the eigen- 
values of gf, and is minimized (maximized) when | ip) = 
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|??mm) {\'nmax))- For any IV') ^ \ rirain) {\'nmax)), there 
exists a state, | ip') , in the neighbourhood of | ip) such 
that v^/ < (v^/ > v^), thus this global optimiza- 
tion problem can easily be solved using standard iterative 
methods, such as steepest decent 

Estimation of extremal eigenvalues plays a significant 
role in the direct detection [Q and distillation [|j of quan- 
tum entanglement. For example, in a special case if two 
qubits described by the density operator g;,, such that 
the reduced density operator of one of the qubits is max- 
imally mixed, we can test for the separability of Qh by 
checking whether the maximal eigenvalue of Qh does not 
exceed j fict] . 

Finally, we may want to estimate an unknown state, 
say a d X d density operator, gi,. Such an operator is de- 
termined by c?^ — 1 real parameters. In order to estimate 
matrix elements ('0 | | "0) i I'un the network as many 
times as possible (limited by the number of copies of gb 
at our disposal) on the input \ip) {ip \ ® gb, where | "0) is 
a pure state of our choice. For a fixed | 0), after several 
runs, we obtain an estimation of, 

V = I Pfc I V) • (6) 

In some chosen basis { | ?^) } the diagonal elements 
(n I gf, I n) can be determined using the input states 
I n) (n I (X) gb- The real part of the off-diagonal ele- 
ment (?T. I I k) can be estimated by choosing | ip) = 
(I n) -I- |fc))/\/2, and the imaginary part by choosing 
I V') = (I ") + * I fc))/ V^. In particular, if we want to esti- 
mate the density operator of a qubit, we can choose the 
pure states, | 0) (spin +z), (| 0) -I- 1 1)) /v^ (spin +x) and 
(|0) + i|l))/V2 (spin +y), i.e. the three components of 
the Bloch vector. 

Needless to say, quantum tomography can be per- 
formed in many other ways, the practicalities of which 
depend on technologies involved. However, it is worth 
stressing that our scheme is based on a network of a fixed 
architecture which is controlled only by input data, a fea- 
ture that can be useful in some quantum communication 
scenarios. 

We can extend the procedure above to cover estima- 
tions of expectation values of arbitrary observables. This 
can be done with the network shown in Fig. because 
estimations of mean values of any observable can al- 
ways be reduced to estimations of a binary two-output 
POVMs [pT| .We shall apply the technique developed in 
Refs. |14| . As A' = ^I + A is positive if —7 is the mini- 
mum negative eigenvalue of A, we can construct the state 

8a — QA' — rp and apply our interference scheme to 
ir (A') 

the pair gA' ^ gb- The visibility gives us the mean value 
of V, 

V = {V)e^,0g, = Tr (^-^^^^gb^ , (7) 
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FIG. 3: A quantum channel A acting on one of the subsystems 
of a bipartite maximally entangled state of the form | ■(/)+) = 
I fc) I fc) /^/d. P+ is the corresponding density operator, 
i.e. P+ — 1 ]. The output is the state described by 
the density operator qa = 2 X)fci \ k) {l \ igi A {\ k) {I |), which 
contains a complete information about the channel. This iso- 
morphism between A and qa allows to infer all properties of 
the channel from the corresponding properties of the state. 
Any subsequent estimations of qa, or any of its functions, 
provides information about the completely positive map A. 

which leads us to the desired value, 

(A),, = Tr {gbA) = vTr^ + 7(vd - 1), (8) 

where Tr I = d. 

Any technique that allows direct estimations of prop- 
erties of quantum states can be also used to estimate cer- 
tain properties of quantum channels. Recall that, from 
a mathematical point of view, a quantum channel is a 
trace preserving linear map, g A(p), which maps den- 
sity operators into density operators, and whose trivial 
extensions, Ife (g) A do the same, i.e. A is a completely 
positive map. In a chosen basis the action of the channel 
on a density operator 9 = '^^i Qki \ k) {I \ can be written 
as 

A{g) = A ( ^ g,i | fc) (/ 1 ) = ^ (| k) {I |) . (9) 

\ kl / kl 

Thus the channel is completely characterized by opera- 
tors A(| fc) (/ 1). In fact, with every channel A we can 
associate a quantum state g\ which provides a complete 
characterization of the channel. For if we prepare a max- 
imally entangled states of two particles described by the 
density operator P+ = i ^^^^^ \k) {l\ \k) and if we 
send only one particle through the channel, as shown in 
Fig. H, then we obtain 

P+ ^ [I A] P+ = gA, (10) 

where 

eA-^Eifc)('i®A(|fc)('i)- (11) 

kl 

We may interpret this as mapping the | k) {I |*^- 
element of an input density matrix to the output matrix, 



4 



A (I k) {I I). Thus, knowledge of qa allows us to deter- 
mine the action of A on an arbitrary state, g A(p). 
If we perform a state tomography on g\ we effectively 
perform a quantum channel tomography. If we choose 
to estimate directly some functions of qa then we gain 
some knowledge about specific properties of the channel 
without performing the full tomography of the channel. 

For example, consider a single qubit channel. Sup- 
pose we are interested in the maximal rate of a reliable 
transmission of qubits per use of the channel, which can 
be quantified as the channel capacity. However, unlike in 
the classical case, quantum channels admit several capac- 
ities jl^, |l^ , because users of quantum channels can also 
exchange classical information. We have then the capac- 
ities Qc where C = 0, stands for zero way, one 
way and two way classical communication. In general, it 
is very difficult to calculate the capacity of a given chan- 
nel. However, our extremal eigenvalue estimation scheme 
provides a simple necessary and sufficient condition for 
a one qubit channel to have non-zero two-way capacity. 
Namely, > iff has maximal eigenvalue greater 
than i. (Clearly, this a necessary condition for the other 
three capacities to be non-zero). 

This result becomes apparent by noticing that if we 
trace gA over the qubit that went through the channel 
A (particle 2 in Fig. |^), we obtain the maximally mixed 
state. Furthermore, the two qubit state, gA, is two-way 
distillable iff the operator | (g) I — has a negative eigen- 
value (see for details), or equivalently, when gA has 
the maximal eigenvalue greater than i. This implies 
Q^{A) > because two-way distillable entanglement, 
which is non-zero iff given state is two way distillable, is 
the lower bound for Q^{A) |T^. 

In summary, we have described a simple quantum net- 
work which can be used as a basic building block for 
direct quantum estimations of both linear and non-linear 
functionals of any density operator g. It provides a direct 
estimation of the overlap of any two unknown quantum 
states ga and gt, i.e. Tr pagb- Its straightforward exten- 
sion can be employed to estimate functionals of any pow- 
ers of density operators. The network has many potential 
applications ranging from purity tests and eigenvalue es- 
timations to direct characterization of some properties of 
quantum channels. 

Finally let us also mention that the controUed-SWAP 
operation is a direct generalization of a Fredkin gate 
and can be constructed out of simple quantum logic 



gates [|. This means that experimental realizations 
of the proposed network are within the reach of quan- 
tum technology that is currently being developed (for an 
overview see, for example, [p^). 
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